Systematic computation of crystal field multiplets for X-ray core spectroscopies 
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We present an approach to computing multiplets for core spectroscopies, whereby the crystal field 
is constructed explicitly from the positions and charges of surrounding atoms. The simplicity of 
the input allows the consideration of crystal fields of any symmetry, and in particular facilitates 
the study of spectroscopic effects arising from low symmetry environments. The interplay between 
polarization directions and crystal field can also be conveniently investigated. The determination 
of the multiplets proceeds from a Dirac density functional atomic calculation, followed by the exact 
diagonalization of the Coulomb, spin-orbit and crystal field interactions for the electrons in the open 
shells. The eigenstates are then used to simulate X-ray Absorption Spectroscopy and Resonant 
Inelastic X-ray Scattering spectra. In examples ranging from high symmetry down to low symmetry 
environment, comparisons with experiments are done with unadjusted model parameters as well as 
with semi-empirically optimized ones. Furthermore, predictions for the RIXS of low-temperature 
MnO and for Dy in a molecular complex are proposed. 

PACS numbers: 78.70.Dm, 78.70.En, 78.70.Ck, 71.70.Ch 



I. INTRODUCTION 

X-ray Absorption Spectroscopy (XAS) and X-ray Res- 
onant Inelastic X-ray Scattering (RIXS) measurements 
act as a local probe for various excitations in materials. 
In particular, the spectra exhibit prominently the fea- 
tures associated to atomic multiplets generated by open d 
and / shells. Brought by the development of synchrotron 
light sources, a lot of progress has been made in the last 
few years in RIXS, a review of which is found in Ref. [l|. It 
is now established that XAS and RIXS, especially in the 
soft X-ray regime^, have become the tools of choice to in- 
vestigate transition metal compounds. Typical examples 
are the titanates 3 - - — , the vanadates^— and the cupratcs 
in the harcU^ and soft X-ray regime o 13 : 14 . Lately the in- 
creased resolving power of the instruments has allowed 
the observation of multiparticle excitations 1 ^— arising 
in inelastic X-ray experiments. 

As the complexity of experiments increases and the 
range of materials investigated expands, it has been felt 
that there is the need for an easy tool relating the energy 
and polarization dependence of the observed spectra to 
the properties of crystal field atomic multiplets. 

It has been known early on that the energy levels of 
the electrons in an open shell form multiplets under the 
effect of the electron-electron interaction, the spin-orbit 
coupling, and the crystal fields. While many aspects of 
the splitting can be discussed from a central field model 
and symmetry considerations alone, a quantitative ap- 
proach from first principles requires the use of computers. 
A review of the pioneering developments in computation 
is given in Cowan's book^i. The early application of 
Cowan's approach concerned the spectra of atoms in the 
ultra-violet and visible range. An important further step 
was the inclusion of relativistic corrections to the atomic 
radial functional, which was used together with a sec- 
ond variation treatment of intermediate coupling involv- 
ing spin-orbit splitting and electron-electron interaction. 



It was realized in the early eighties that atomic core 
spectroscopies such as X-ray photoemission spectroscopy, 
XAS, RIXS and similar spectra can show multiplct struc- 
ture signatures when open d- or f- shells are involved^ 3 . - — . 
An early theoretical study of atomic multiplct splitting in 
a cubic crystal field has been published by van der Laan 
et a@£, followed up by de Groot et a$L for 3d transition- 
metal compounds. 

Spectroscopic studies are typically applied to solid 
state compounds and therefore should in principle take 
into account hybridization effects. Parameterized model 
Hamiltonians are used to couple levels with a conduction 
band^S and to include charge transfer from the selected 
configurations of the ligands^. The inclusion of charge 
transfer effects in multiplet calculations is required in or- 
der to reproduce with accuracy the spectra of compounds 
with high ligand hole character—. It has however been 
pointed out that the simple crystal field model explains 
a large range of experiments 3 -^. The present study will 
concentrate solely on the effects of the crystal field of any 
symmetry. 

Atomic multiplct codes based on Cowan's pioneering 
work 3 ^— make intense use of symmetry considerations. 
Starting from the LS atomic symmetry standpoint, the 
intermediate coupling arising from spin-orbit coupling is 
included, followed by the point group symmetries for the 
crystal field. The crystal field is entered in parameterized 
form, which depends on the point group in consideration. 
A somewhat simpler and more general method, where the 
crystal field is specified via parameters in a matrix, has 
also been reported 3 ^. Recently, Zhang et rationalized 
counter-intuitive crystal field fit parameters using a point 
charge model with variable ion positions. 

In the present paper we focus on intra atomic electron- 
electron interactions of relativistic atoms in a general 
crystal field and its effect on XAS and RIXS spectra. The 
crystal environment of the atom of interest, rather than 
being characterized by symmetrized matrix elements, is 
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entered explicitly via the positions and charges of the 
neighboring atoms. The input is therefore straightfor- 
ward and requires no prior knowledge of point group sym- 
metry. The multiplet levels are calculated by determin- 
ing the electron-electron and spin-orbit interactions, as 
well as the effect of the crystal field point charges for the 
relevant core and valence shells in a determinantal basis. 
The spin-orbit matrix elements and the radial wave func- 
tions are obtained from the solution of the atomic central 
field Dirac equation in the density functional theory for- 
malism. The resulting eigenvalue problem is solved using 
standard diagonalization methods. 

The paper is organized as follows. The approach to 
calculating multiplets with the point charge model defin- 
ing the crystal field is presented in Sec. [TTJ The XAS and 
RIXS spectra most commonly studied today are domi- 
nated by dipolar transitions caused by the interaction of 
polarized light with the relevant core and valence elec- 
trons of the atom. A brief reminder of the standard the- 
oretical treatment for XAS and RIXS and the selection 
rules are given is Sec. IIIII The application of the present 
approach to both well known high symmetry and less fa- 
miliar low symmetry cases will be demonstrated is Sec. 
|TV]for XAS and in Sec. |V]for RIXS. Summary, conclu- 
sions and outlook follow in Sec. IVTl 



II. 



ATOMIC MULTIPLETS FOR GENERAL 
CRYSTAL FIELD 



The general theory of atomic multiplet is described ex- 
tensively in Refs. [2(| Hi], Hi], and H3 Within the point 
charges model, the relevant interactions for the electrons 
in the open shells are the intra-atomic electrostatic in- 
teractions, the spin-orbit coupling and the effect of the 
crystal field. Considering n electrons in the open shells 
of an atom and N- lons point charges Q m at position R r „ 
from the atom, the multiplet Hamiltonian is given by 



electrons in 2p and in the 3d shells will be included in 
Eq. HJ 

The bulk of this part of the paper addresses some as- 
pects of the implementation of the three terms of Eq. 
[1] The Hilbert space is presented first. Then, starting 
with the central field atomic equation, we discuss the 
spin-orbit coupling, the crystal field, and the effect of 
ionization. A comparison with other existing approaches 
is conducted at the end the section. 



A. The Hilbert space 

The focus of this paper lies in calculating multiplets 
involved in soft X-ray absorption and emission processes. 
The wave functions are antisymmctrized products of the 
spin-orbitals of the electrons in the atomic open shells, 
so the dimension of the Hilbert space is, in the present 
model, 



AT(l, k) 



[2(2/ + 1)]! 



[2(2Z + !)-&]! k\ 



(3) 



for a shell labeled by the quantum numbers (n, /) and 
containing k electrons. For instance 5 electrons in a p- 
shell gives A/"(l,5)=6 determinants, 2 electrons in a d- 
shell gives Af(2, 2) = 45 determinants. If more than one 
shell is opened, each shell being independent, the overall 
size is given by the product Af = Yl^i- The case of 
transition metals can be handled easily on today's stan- 
dard desktop computers. The largest Hilbert spaces in 
present day soft X-ray applications are found in the case 
of lanthanides, for which d — > f transitions with an open 
/ shell occur. An example could be Gd 3+ in the ground 
state configuration 3d 9 4/ 7 , with a final state matrix di- 
mension of 34320. The generation and diagonalization 
of this size of matrices remain however still tractable on 
contemporary architecture with sufficient memory. 
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The third term V x t a i is the crystal field interaction in the 
crystal field potential V x t&i 
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The first term in Eq. [T]is the electron-electron Coulomb 
interaction. The second is a single electron operator ac- 
counting for the kinetic energy and spin-orbit interaction; 
this term will be drawn from the solution of the atomic, 
central field Dirac equation, a point that will be discussed 
further below in this section. A typical case involving a 
transition metal will require considering the multiplets 
arising from opening a hole in the 2p shell and promot- 
ing an electron in an incomplete 3d shell, so that all the 



B. Relativistic electronic orbitals and energies for 
a single atom 

Computing atomic multiplet levels requires central 
field model orbitals and energies. Previous atomic multi- 
plet calculations, including the more recent approaches^, 
are usually based on the non-relativistic Schrodinger 
equation. In this work the Dirac equation has been pre- 
ferred, since it naturally includes spin-orbit coupling ef- 
fects. 

In the absence of crystal field, the Dirac equation is 
written in a formal way as follows: 

H D * ! = e^i ; U D = cot ■ p + (3mc 2 + V(r) (4) 



a = 



er 
<t 



h 






-]1 2 



(5) 



3 



where <r are the Pauli matrices and V(r) is a spherically 
symmetric field which can be decomposed as follows: 

V(r) = -- + V s {r) + V xc (r) ; (6) 
r 

The first term is the electrostatic attraction by the nu- 
clear charge Z and Vg is the static potential: 

r oo 

V s (r) = - [ p(x)4nx 2 dx + [ ^4iTx 2 dx. (7) 
r J J x 

r 

In the multiplet code developed for the present work 
the Dirac equation is solved for the neutral atom (see 
Sec. Ill C|) . within the density functional theory^ and us- 
ing a local density functional approximation V xc {r). The 
radial wavefunctions and eigenvalues Si will be used to 
compute the matrix elements of Eq. [TJ For light ele- 
ments, the non-relativistic limit of the radial wave func- 
tions can be retrieved by performing a weighted average 
of the radial term over the shell. For the examples shown 
in this work the crystal field matrix elements have also 
been calculated in this limit. 

For the application to XAS and PJXS the Dirac equa- 
tion has to be solved twice: once in the ground state 
and once with an electron moved from a core shell to a 
valence shell. 



C. Spin-orbit coupling and Coulomb term 

As mentioned above, the eigenvalues Ei of the Dirac 
equation are used to build a diagonal operator that, by 
default, accounts for the effect of spin-orbit coupling (sec- 
ond term in Eq. QJ). All the electrons of the participating 
core and valence shells contribute to the spin-orbit oper- 
ator. 

Such an operator includes the kinetic energy, but also 
contains the contribution of the central symmetric part 
of the electron-electron interaction, which in principle is 
already addressed by the first term (Coulomb term) in 
Eq. [TJ 

Double counting can be partially prevented by consid- 
ering the following standard expansion that applies to 
the Coulomb term in Eq. [TJ 

1 k—0 m— — k 1 

oo (8) 

with n < T2- Replacing the sum X)fc=o w ith 2fc=i m 
the expansion of Eq. [5] effectively removes a spherically 
symmetric term corresponding to V s in Eq. and part 
of the exchange interaction V xc . 

Alternatively, and this is the approach used in all the 
examples presented in this paper, the radial term in the 
matrix element of the Coulomb interaction is averaged 
over the shell, thereby effectively removing unwanted 
contributions. 



TABLE I. Positions and charges of the crystal field for Cu 2+ 
in monoclinic tenorite CuO. The oxygen coordinates x,y,z and 
distance r to Cu are given in A, the charges in units of e. For 
comparison, the cubic polymorph of CuO has 6 oxygen atoms 
at a distance of r=2.12 A. 



x 


y 


z 


r 


charge 


1.3304 


0.5617 


-1.2461 


1.910 


-2.0 


0.9001 


1.1298 


1.2836 


1.930 


-2.0 


-1.0143 


-1.1517 


-1.2461 


1.980 


-2.0 


-1.4446 


-0.5836 


1.2836 


2.020 


-2.0 


0.9001 


-2.2970 


1.2836 


2.780 


-2.0 


-1.0143 


2.2751 


-1.2461 


2.790 


-2.0 



It will be seen from comparisons with X-ray spectro- 
scopic experiments in Sec. IIVI and [V] that this straight- 
forward approach to spin-orbit coupling delivers a good 
first-principle prediction of the L3-L2 edge splitting. 



D. Crystal Field as point charge ligands 

Unlike group-theoretical approaches to calculating 
multiplet levels, the present work accounts for crystal 
field splitting by explicitly including the electrostatic in- 
teractions of neighboring point charges with the open 
shells of the atom. The crystal field is then computed 
together with the direct and exchange Coulomb interac- 
tions, and the spin-orbit coupling (Eq. [TJ). 

The advantage of treating the crystal field as point 
charge ligands is two-fold: i) the position of each ion is 
given as input in the program without having to resort to 
symmetry and ii) the crystal field arising from charges at 
experimental atomic positions leads to a reasonable pre- 
diction with often sufficient accuracy for interpretation 
without using a fitting parameter. 

The purpose of the charges is to recreate the proper sym- 
metry of the crystal. It is natural to take as a starting 
point the oxidation state, or the formal charges, the lat- 
ter reflecting the level of band filling in insulators. Usu- 
ally however, either fail to reproduce correctly the actual 
charge distribution in the solid, and it would be legiti- 
mate to treat the charges assigned to equivalent atoms 
as fitting parameters. This exercise has not been carried 
through in the examples discussed in this work. 

An example of input is given in Table [I] for the mono- 
clinic tenorite (CuO), where the Cu has 6 nearest neigh- 
bors forming a distorted octahedral environment. The 
positions of the neighbor O atoms are taken from the In- 
organic Crystal Structures Database (ICSD) entry num- 
ber 69758. A formal charge of -2e was allocated to each O 
atom. The ground state configuration of Cu 2+ is 2p 6 3d 9 . 
This information, together with the elements of Table [TJ 
are the only input required in the code. In that par- 
ticular case, and in the non-relativistic (NR) limit, the 
splitting of levels originates solely from the crystal field. 
The output of the code is represented graphically in Fig. 
[JJ The first column on the left shows the levels for Cu 2+ 
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(atomic units) 



FIG. 1. Multiplet levels of Cu 2+ in CuO calculated with the 
present code in the NR limit. Column 1: cubic octahedral 
environment (6 O atoms). Columns 2-5: monoclinically dis- 
torted environment built up of the 6, 16, ~ 500 and ~ 1600 
closest ions, respectively. 

in octahedral symmetry. This structure exists and is a 
polymorph of CuO in a cubic phase (ICSD 61323), where 
the six O atoms are forming an octahedron at a distance 
of 2.12 A to the Cu. If the crystal field of Tablets con- 
sidered instead (column 2), the energy levels originally 
ascribed to the tig (3-fold degenerate) and e g (2-fold de- 
generate) orbitals in the cubic case split into three, re- 
spectively two further levels. 

The convergence with increasing numbers of ions as lig- 
ands is illustrated by the columns 2-5 in Fig. [TJ with 
respectively 6, 16, s» 500 and ~ 1600 closest ions consid- 
ered. While the explicit symmetry breaking induced by 
the crystal field is already accounted for by the 6 closest 
oxygen atoms, the addition of the second shell of ions 
(16 closest ions) is necessary to reach a better qualita- 
tive agreement with the limit. In the particular example 
of monoclinic CuO, the convergence is achieved for most 
practical purposes with about 500 ions. 

E. Semi-empirical parameters 

The assumption behind the crystal field model is that 
the atom under consideration is surrounded by point 
charges and is itself in an particular ionized state. The 
ionization has a major impact on the radial part of the 
wave functions of an atom. Moreover, the radial wave 
functions will be affected by the bonding environment of 
the atom. The latter consideration is the origin of the 
introduction of semi-empirical parameters. As the radial 
functions enter in both the calculation of the crystal field 
and the Coulomb term, the sensitivity of these terms to 
the formal ionization state is addressed below. 

An atom give rises to different radial wave func- 
tions according to the degree of ionicity. The reduced 
screening felt by the remaining electrons of a positively 



FIG. 2. (Color online) Radial part of the d-wave function for 
Ti as a function of r in atomic units: the tail of the wave 
function shrinks with increasing ionization. 



charged atom leads to the contraction of the radial 
wave functions around the nucleus. As illustrated in 
Fig. [2] for the case of Ti, the expectation value (r) 
in the Ti 3+ ion varies by about 30% compared to the 
neutral titanium ion. An estimation of the effect of 
the contraction of the radial function can be given as 
follows. For an open d shell in a non-relativistic limit, 
the only contribution to level splitting comes from the 
anisotropic densities with a purely k = 4 multipolar 
character. Eq. [5] implies that the crystal field splitting 
scales with the power of the anisotropic density. As the 
radial wave function becomes sharper with increasing 
ionization, the crystal field splitting for a shell of angular 
momentum I scales as (r +21 ). On the other hand, the 
electron-electron interaction scales like (r _1 ). As a 
consequence, if the wave function does actually extend 
further than anticipated, the crystal field strength would 
be under-evaluated, whereas the Coulomb interaction 
would be over-evaluated. 

A comparison with experimental spectra indicates that 
a better correspondence is achieved if the radial func- 
tions of the neutral atom with integral ionic charge num- 
ber is used, rather than those of the positively charged 
ion. Moreover, self-consistent band-structure calcula- 
tions suggest that the charge density around atoms in 
a solid is close to that of neutral atoms. All the exam- 
ples presented in this paper have therefore been obtained 
with radial functions calculated from the neutral atom. 

Further comparisons with experiments show, unsur- 
prisingly, that an even better match can be made if 
the matrix elements of the Coulomb interaction, the 
spin-orbit coupling and the crystal field are allowed 
to vary. A scaling of these elements can be seen as a 
way to adjust the radial functions to the hybridization 
environment encountered in the solid. As will be seen 
in examples below, it is usually necessary to scale the 
crystal field up, and the coulomb interaction down. From 
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the discussion in the paragraph above, this suggests that 
the neutral atom radial wave functions still fall off too 
rapidly. 

An empirical correction to the ionic model is therefore 
done in our code via the parameters S cou i, S x tai and S soc , 
which scale the respective terms in the Hamiltonian 
(Eq. [T]). The crystal field scaling may also incorporate a 
global adjustment to the chosen charges. 

We will call "raw" the multiplet results or spectra 
calculated using formal charges, respectively oxidation 
states, and crystallographic neighbor positions, all 
scaling parameters being set to unity. 

The model of a single ion in an ionic environment ig- 
nores charge transfer from ligands, band effects and col- 
lective excitations. It has nonetheless been successfully 
applied to the calculation of multiplct-dcpcndent spec- 
tral features 3 ^. Ghiringhelli et al compared the crystal 
field model to the single impurity Anderson model for a 
prototype transition metal oxide, MnO, in a cubic en- 
vironment. They found that both models agree well in 
that part of the spectra determined by local effects. The 
multiplet approache a 26 i 27 based on Cowan's workSi use 
Hartrcc-Fock parameters (which include relativistic cor- 
rections) . These values are usually scaled down by 70 to 
80% to compensate for the neglect of configuration inter- 
actions in the Hartree-Fock theory. The effective crys- 
tal field parameters are adjusted to experiments. The 
present approach makes use of the DFT-LDA results of 
the relativistic Dirac equation; in principle it delivers the 
multiplet structure without the need of parameterization; 
if fitting to experiments is required, the scaling parame- 
ters and possibly the charges can be adjusted. 



F. Results and comparison 

In the presence of a cubic crystal field and in the non- 
relativistic limit, the energy levels of one electron in a d 
atomic orbital split into a four-fold degenerate e g and a 
six- fold degenerate t 2g orbital (including spins). In this 
case, the cubic crystal field splitting is well defined as 
the energy difference between the e g and the t 2g lev- 
els. This value is also often named 10D q in the liter- 
ature, and is taken as an effective parameter fitted to 
experiments 2 ^. In Ref. H3, the authors estimated that for 
transition metals, the crystal field splitting would typi- 
cally vary between and 2.5 cV. They point out that the 
splitting found by examining optical spectroscopy exper- 
iments correspond to the final state splitting, which can 
in principle differ from that of the ground state initial 
value. 

As the present work deals with low as well as high 
symmetry, a characteristic crystal field splitting (CXS) is 
defined, which expresses the strength of the crystal field 
splitting. For cubic symmetry the CXS is equivalent to 



TABLE II. Multiplets of the ground state 2p 6 3d 1 of LaTi0 3 . 
Eigenvalues are given in eV for unsealed (raw) contributions 
to the Hamiltonian (Eq. [T]), and with scaling determined 
from fitting to XAS data 4 ^, as explained in the text. The 
degeneracy deg (including spin) is indicated for the cubio-S 
and orthorhombia^i ligands field. 





Cubic 
deg raw 


Orthorhombic 
deg raw scaled 


Crystal field 
(NR limit) 


6 -1.083 
4 1.624 


2 -1.200 -1.500 
2 -0.849 -1.062 
2 -0.735 -0.919 
2 1.327 1.659 
2 1.458 1.822 


Spin-orbit 


4 -0.027 
6 0.018 


4 -0.027 -0.025 
6 0.018 0.017 


Crystal field 
+ Spin-orbit 


4 -1.092 
2 -1.065 
4 1.624 


2 -1.200 -1.500 
2 -0.850 -1.062 
2 -0.735 -0.919 
2 1.327 1.659 
2 1.458 1.822 



the cubic crystal field splitting. In the present approach, 
the CXS can be easily determined for a cubic field by 
evaluating the central atom in a d 1 or d 9 configuration 
in the non-relativistic limit. In the example of Cu in 
cubic geometry, as can be seen in Fig. [TJ the unsealed, 
raw result of the calculation gives CXS=0.6 eV. In the 
monoclinic structure the strength of the crystal field can 
be taken as the energy difference between the barycenter 
of the two lowest levels and that of the 3 highest levels, 
which gives 1.9 eV. Alternatively, one can take the CXS 
as the difference between the highest of the two lowest 
and the average of three highest levels, which amounts 
to 0.7 eV (for the « 500 ions case). 

A well studied example is Ti 3+ in LaTiOs, in its 2p e 3d 1 
ground state configuration. Once thought of as a perfect 
cubic perovskite 4 ^, LaTiOs has been characterized a 
number of times; a later powder diffraction study 41 
detected small orthorhombic distortions to the octahe- 
dral environment. The multiplet results of our code for 
LaTi0 3 in the idealized cubic structure^ (ICSD 28908) 
and with orthorhombic distortions 41 - (ICSD 98414) is 
presented in Table [TTJ The formal charges La 3+ , Ti 3+ 
and O 2- were allocated to the ions. The raw, "first- 
principles" eigenvalues are obtained with, as sole input, 
the specification of Ti as the atom under consideration, 
its ground state configuration 2p 6 3d 1 and the crystal 
field as described above. For both structures all charges 
within 20 A of the Ti atom were included. The first row 
is obtained by considering only the Coulomb interactions 
and the crystal field in the NR limit. As expected, the 
cubic crystal field causes a t2 g -e g split, giving CXS=2.7 
eV. The effect of the spin-orbit coupling, without the 
crystal field (second row in Table HI]) , is to split levels 
by 45 meV. The combination of cubic crystal field and 
spin-orbit coupling has the effect of splitting the t 2g 
level by 27 meV (third row in the Table). With the 
orthorhombic field (column labeled Orthorhombic raw) 
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all the levels are already split by the crystal field in 
the NR case. The inclusion of spin-orbit interactions 
(third row in the figure) has only a minute effect on the 
levels. Taking CXS as the energy difference between 
the highest negative value and the average of the posi- 
tive value, we find CXS=2.1 eV in the orthorhombic case. 

The scaled results for LaTiC>3 are shown in Tabic [XT] 
along the original, raw values for the orthorhombic ge- 
ometry. The scaling parameters were obtained by fitting 
to XAS experiments^ (see results in Sec. IIV|) . which 
gave S cou i=0.7, S x tai=T-25 and S SOC =0.915. The cor- 
rected CXS now has a value of 2.7 eV (2.1 eV unsealed). 
In the scaled case the splitting of the ti g level is of about 
0.14 and 0.4 eV. These values agree well with the mea- 
surements and LDA-based, Wannier function projection 
calculations of Ref. [43. 



III. XAS AND RIXS THEORY 

XAS and RIXS involve processes whereby an X-ray 
photon is absorbed by a core electron, which is promoted 
into a valence shell. The resulting core-hole and valence 
configuration, in general, gives rise to a multiplet struc- 
ture. In the present work the ground state and core- 
hole state are determined using the multiplet Hamilto- 
nian (Eq. [1]) as sketched in Sec. HU 

XAS measures the decay products of the core-holei 
In RIXS measurements the system returns to the ground 
state by emitting a photon. The standard theoretical 
formulation for the calculation of the signal intensity of 
XAS and RIXS is stated below, as well as the dipolar ap- 
proximation used in both cases. Lastly, a brief reminder 
of the polarization effects and selection rules is given. 



A. XAS intensity formula 

The XAS absorption spectra is simulated the usual 
way^I with the Fermi Golden rule 



I(w) cx Y,\(^P\^o)\ 2 S{hw + E - E t ) 



(9) 



fajj is the energy of the absorbed photon. |?/>o) refers to 
the ground-state and Eq to its eigenvalue, which is often 
degenerate. \tpi) and Ei are respectively all the possible 
eigenvectors and energies of the excited, core-hole state. 
The transition operator O is treated in the dipolar ap- 
proximation, briefly exposed in Sec. IIII C[ and depends 
on the polarization of the incoming light e,- n (Sec. IIII D[) . 

Intrinsic effects as well as the finite experimental res- 
olution leads inevitably to the broadening of the sharp 
delta peak of Eq. |Hl The spectral broadening here is 
approximated by a Lorentzian. Eq. [5] becomes: 



An intrinsic contribution to the Lorentzian broadening 
comes from the finite lifetime of the core-hole state. 
A typical £2,3 core-hole linewidth for transition metals 
varies between 0.2 and 0.4 eVi£. The different level of 
hybridizations of the states involved may however be re- 
flected in a variation of the vibrational and dispersional 
broadening of the peaks4^. 

If F; = F is taken as a constant, the relative intensities 
of the simulated peaks may differ from that observed in 
experiments. A Gaussian experimental broadening can 
be applied as well if a better fit to experiment is sought; 
this has not been done in the examples shown in this 
work. 

The multiplet levels of the excited core-hole state are 
given with respect to the energy difference between the 
total energy of the single neutral atom in the ground 
state and in the excited state, which in the present case 
is determined within the LDA-DFT approximation. The 
absolute energy position of the multiplets is not deter- 
mined with sufficient accuracy and, compared to experi- 
ments, tends to shift the peaks to a lower energy by an 
amount of the order of the percent. It is usual to align 
the calculated spectrum to the measured one, and this 
has been systematically done in this work for the XAS 
spectra. 



B. RIXS intensity formula 

RIXS processes involve the creation and subsequent 
radiative deexcitation of a core hole. With the incoming 
photon energy huJi n close to an absorption edge, the RIXS 
intensity is approximated, using the Kramers-Heisenberg 
formula^, by the resonant term 



E 



v ^ f \O A \^i){^\O\i> ) 
^ E,, 



Eq — fkOin — iTi 

xS(Ef-Eo-nn) (11) 



I(u>) (x]T|<^|O|V>0>| 5 



{hw + Eo-Erf + r? 



(10) 



Details on the second order perturbation theory leading 
to Eq. [TT]can be found in the review Ref. QJ. h£l is the en- 
ergy loss, i.e. the energy difference between the incoming 
and outgoing photon. During the first stage of the pro- 
cess a core-electron is promoted from the ground-state 
|^o) °f energy Eq into an excited, core-hole intermediate 
state \ipi) of energy Ei. The probability of this transition 
depends on the polarization of the incoming light and is 
determined by the dipolar operator O. The deexcitation 
from the intermediate to a final state \ipf) of energy Ef 
occurs via the emission of a photon of energy hujf and 
depends on the polarization of the outgoing photon, as 
given by the operator O' (see Sec. IIII D[) . The core- hole 
lifetime Ti, as above in the XAS case, is often taken as a 
constant. The S function in Eq. [11] is also broadened by 
the finite lifetime of the final state \ipf) and experimental 
resolution. A Lorentzian function is usually used in lieu 
of the delta function, as is the case for the present work. 
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The experimental RIXS spectra also contain broad flu- 
orescence lines which mainly come from the emission 
from intermediate states with finite lifetimes which are 
not taken into account in the multiplet calculation. Non- 
local contributions are also ignored in this study for the 
present. 

C. Dipolar Approximation 

The incoming photon is defined by its polarization e 
and its wavevector k. After expanding the vector poten- 
tial in plane waves, and in the limit k • r <C 1, the matrix 
elements in Eq. [5] and [TT] can be obtained as 

(^:|e-pe lk - r |V'o) ~ (ipi\e • P (1 + *k • r - ...) |Vo) (12) 

The dipole approximation consists in keeping only the 
first term, such that the operator is given by 

C> = e-p = e- [r,H] oc e • r (13) 

Using an expansion on the F™, we finally obtain 

6 cx r (e^ 1 + e Y? + t-iYf 1 ) (14) 

where the coefficients Ci represent the projection of the 
polarization vectors on the Y™ basis. 

At the dipolar approximation level the selection rules 
for the transition between states are the following: 

A/ = ±1 ; As = ; A J = 0, ±1; ... (15) 

In the present code these rules are automatically included 
through the computation of Gaunt coefficients formed by 
combining the spherical harmonics of Eq. I14l and that of 
the wave functions of the matrix elements in Eq. [9] and 
ITT1 The weights of the Gaunt coefficients are given by the 
polarization projections. The polarization has therefore 
a discriminating effect on the transitions. 

D. Polarization effects 

The approach presented in this work is well suited to 
the study of polarization effects in low symmetry. The 
polarization of the photon is given in the same frame of 
reference as the ions forming the crystal field, and there- 
fore any possible geometries can be analyzed quantita- 
tively. An example is shown in Fig. [3] for an atom in 
a planar environment. The 4 ions forming the crystal 
field are expressed in the local (x, y, z) basis, and so is 
the polarization of the incoming and outgoing photon. In 
most RIXS experiments the polarization of the outgoing 
photons is currently not recorded, and only the detec- 
tion direction £> out of the beam is known. The signal is 
simulated by the incoherent superposition of intensities 
contributed by the polarizations orthogonal to the detec- 
tion direction. A similar superposition is done in the case 




X 



FIG. 3. (Color online) Polarization directions with respect to 
the crystal field: an incoming beam of photons of polarization 
ein excites the central atom (small purple ball). The outgoing 
photons that are detected along a direction 6 out can have any 
polarizations perpendicular to it. As an example here ej} ut and 
e^ ut are also parallel, respectively perpendicular, to e ln . The 
four ions (large red balls) making up the planar environment 
of the central atom are described in the same reference frame 
(x, y, z) as the polarization and beam directions. 

of unpolarized photons, or experiments done on powder 
samples. 

As pointed out above, the weights of the Gaunt coef- 
ficients are determined by the polarization. In the most 
general case different polarizations will suppress or en- 
hance other transitions within the dipole approximation. 
Examples of the effect of polarization directions will be 
given in the Sec. IIVI and fVl 



IV. XAS SPECTRA 

In this section several XAS spectra simulated with 
the present code are compared with experiments. First 
the familiar high symmetry case of the SrTi03 is 
reproduced. Then the spectrum of LaTi03 with cubic 
and with orthorhombic distortions is presented. Finally, 
the monoclinic epidote is discussed. 



A. Cubic SrTiC-3 

SrTiOs is an example of cubic pcrovskitc with formal 
charges Ti 4+ , Sr 2 + and Ti 2 ~. The XAS spectrum at 
the titanium L-edge will therefore involve the transitions 
2p 6 3d° to 2p 5 3d 1 . While no multiplet structures appear 
in the ground state, the Hilbert space for the core-hole 
state involves 60 determinants. The multiplet levels re- 
flects the non-trivial interplay between the intra-atomic 
Coulomb interaction and spin-orbit, and the crystal field, 
as is illustrated in Fig. EJa), a raw result of the code. An 
artificially small core-hole width of 0.05 has been cho- 
sen in order to exhibit the positions of the peak more 
clearly. The cubic crystal field is such that the 6 oxy- 
gen neighbors are at a distance of d = 1.95 A. In the 
non-relativistic limit the pure Coulomb peak (curve 1, 
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FIG. 4. (Color online) XAS spectra for Ti 4+ in SrTiOa in the 
2p 6 3d° ground state configuration, (a) Raw results of the mul- 
tiplet code. Lines from top to bottom: Coulomb interaction 
in the NR limit (1, brown), then adding crystal field (2, pur- 
ple), then for Coulomb interactions and spin-orbit coupling (3, 
green), and finally for Coulomb, spin-orbit and crystal field 
(4, blue), (b) Top line (blue): raw results of the multiplct 
code with all contributions and a constant core-hole lifetime; 
middle line (black): experimental result of J. Schlappa et aZ.— ; 
lower curves (red): results with rescaling of all contributions 
as explained in the text. The dashed line is produced with a 
constant core-hole broadening, the full line with a broadening 
increasing with increasing energy. 



brown) splits in two (curve 2, purple) upon the addi- 
tion of the crystal field. In the absence of crystal field, 
the spin-orbit coupling splits the peak in three (curve 3, 
green), as can be seen with the help of the zoom in the 
inset in the figure. The bottom line (curve 4, blue) is 
the result after including all contributions. Fig. H^a) has 
been produced from the raw results of the code without 
applying any rescaling. Fig. EJb) compares the XAS in- 
tensity calculations with the experimental measurements 
(middle curve, black) taken from Schlappa et a&. This 
measurement is in very good agreement with that of other 
authors^. The top curve (blue) is the raw result of the 
code; it is the same as the bottom curve of (a), but with 
a constant Ti broadening of 0.2 cV, corresponding to a 
realistic intrinsic core-hole lifetime value. The bottom 
curves (red) are the result of rescaling with S cou i=0.85, 
S x tai=l-3 and S SOC =0.93. Keeping the core-hole broad- 
ening constant at 0.2 eV produces the dashed line. In 
this particular case the rescaling of the spin-orbit inter- 
action has the effect of bringing the L3 and L2 peaks 
slightly closer together. The relative heights of the peaks 
is sensitive to the scaling of the intra-atomic Coulomb 
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FIG. 5. (Color online) XAS spectra for Ti 3+ in LaTiOs in the 
2p 6 3d 1 ground state. Top (blue) lines: raw multiplet results 
for cubic crystal field (dashed), orthorhombic crystal field 
with polarization A (full) and B (dotted line) as described 
in the text. Middle (black) line: experiment as reported in 
Haverkort et ai— . Lower (red) lines: rescaled multiplet results 
for the cubic crystal field (dashed), and for the orthorhombic 
crystal field with polarization A (full) and B (dotted line). 

interaction. The energy difference between the two L3 
peaks (and the two L3 peaks) depends predominantly on 
the crystal field strength and therefore determines the 
rescaling of the crystal field. The agreement with exper- 
iment can be further improved by modeling the energy 
dependence of Ti. The inclusion of a theoretical approach 
for the core-hole lifetime is beyond the scope of this pa- 
per, but the broadening can be empirically matched to 
the experiment. The full bottom line (red) in Fig. |U[b) 
was obtained by increasing T,; from 0.2 eV to 0.8 eV lin- 
early with Ei between 457 and 465 eV. 

Our results compare very well with simulated spec- 
tra produced by multiplet approaches where the crystal 
field is introduced by symmetry group parameters and 
strength^. 



B. Cubic and orthorhombic LaTiC>3 

The ground state multiplcts for the Mott-Hubbard in- 
sulator LaTiO^ has already been discussed in Sec. Ill Fl 
The resulting XAS spectra are plotted in Fig. [5] The 
XAS spectra taken by different authors^^ present some 
variations, possibly reflecting the presence of impurities 
in this compounds. The middle curve in Fig. [5] show 
the latest work^. Both the crystal field corresponding 
to the cubic (dashed line in the figure) and to the or- 
thorhombic structure were tested. In the latter case the 
spectrum becomes sensitive to the polarization direction. 
Two different polarizations were chosen for illustration 
purposes. The first one, labeled A (full lines in Fig. [5]), 
coincides with a Ti-Ti direction for a Ti belonging to 
the next shell of Ti atoms. The second polarization di- 
rection, labeled B (dotted lines in the figure), is along 
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FIG. 6. (Color online) XAS spectra for Fe 3+ in monoclinic 
epidote in the ground state 2p 6 3d 5 . The dotted and full lines 
refer to respectively the xy and z polarization. Top (blue) 
lines: raw multiplet code results. Middle (black) lines: ex- 
perimental results of Henderson et aJ-£. Bottom lines: code 
results with rescaled values as given in the text. 



line to the z polarization. The middle curves in the 
figure are the measurements of Henderson et the 
top curves the raw output of the program and the 
bottom the results of rescaling. In this case the following 
scaling parameters were used: S CO ui=0.7, S x t a i=l-8 and 
S SOC =0.915. For the bottom curves the broadening was 
increased from 0.55 eV to 1.2 eV from 708 eV onwards. 

From the three cases presented, one can see that the 
general aspect and relative trend exhibited by the peaks 
for different polarizations are well reproduced by the raw 
output of the code. For the compounds studied here a 
good agreement with the experiment could be achieved 
with the screening of the Coulomb interaction, a slight 
adjustment of the spin-orbit coupling strength, and an 
enhancement of the crystal field strength. The rescaling 
of the crystal field is in effect an adjustment of the ra- 
dial functions, i.e. their extent through the solid. The 
present scaling is therefore consistent with the modifica- 
tions required to bring the point charge model closer to 
a more realistic crystal model. 



the crystallographic b direction of the compound. These 
two polarizations were chosen for their contrasting be- 
haviors. The top curves in Fig. [5] show the raw results 
of the multiplet code, the bottom after rescaling with 
S C oui=0.7, S xta i=l-25 and S soc — 0.915. The raw result is 
calculated assuming a constant core-hole broadening of 
0.25 eV. The scaled results are obtained with a constant 
core-hole broadening of 0.5 eV. The fits were done using 
the orthorhombic crystal field. Fitting starting from the 
cubic crystal field turned out more difficult in view of the 
positioning of the middle peak; the result of the rescal- 
ing using the parameters stated above but applied on the 
cubic crystal field are also shown in the figure (bottom 
dashed line). 



V. RIXS SPECTRA 

Very few high resolution RIXS spectra have been 
published to this date and even fewer relate to structures 
of very low symmetry. We have not been able to find 
a RIXS measurements where the effect of the low sym- 
metry of the crystal field could be made unquestionably 
clear. We therefore demonstrate the RIXS capabilities 
of our code on two high symmetry examples for which 
experimental comparison is available, MnO and NiO. 
Furthermore, the case of distorted MnO is addressed, 
and results for a lanthanide compound are presented. 



C. Monoclinic epidote 

The monoclinic epidote demonstrates the applica- 
bility of the code to low symmetry. Epidote has the 
generic formula Ca 2 Al 2 (Fe,Al)(Si04)(Si 2 07)0(OH). 
The Fe 3+ in this compound is in a distorted octahedral 
environment formed by 6 oxygen atoms. XAS spectra 
were taken by Henderson et a&Z- on an epidote of 
unspecified composition for two polarizations, one in 
the approximate square-plane of the oxygen atoms and 
the other parallel to the axial compression, named xy 
and z polarization respectively. Many epidote structural 
data have been published for different Fe to Al ratio, 
but few differ substantially to one another when used 
as crystal field in the program. The latest neutron 
diffraction measurements^ was chosen for the results 
presented here (ICSD 168464)i2. The XAS spectra for 
the transitions 2p 6 3d 5 to 2p 5 3d 6 are plotted in Fig. 
The dotted lines refer to xy polarization, the full 



A. Cubic MnO 

The Mn L-edge measurement in cubic MnO is a pop- 
ular case for testing the resolution of the RIXS instru- 
ment. The valence state of Mn is Mn 2+ and corre- 
sponds to the ground state 2p 6 3d 5 . The RIXS process 
2p 6 3d 5 -> 2p 5 3d 6 (2p 6 3d 5 )* offers therefore insight 
into a non-trivial and rich multiplet structure. In a first 
place we assume that MnO has the cubic, NaCl structure, 
with a Mn-0 nearest neighbor distance of 2.22 A. A se- 
ries of measurements for several different incoming pho- 
ton energies and two different polarizations were taken by 
Ghiringhclli et al^. Fig. [7] (a) and (b) reproduces two 
sets of results, at 0.7 eV and 3.7 eV above the L3 edge 
(label C and E in Ref. l39h respectively. In this figure, the 
curves measured or calculated with the incoming photon 
polarization horizontal, respectively vertical to the scat- 
tering plane are represented by dotted lines, respectively 
full lines. We refer the reader to Ref. for the descrip- 
tion of the geometrical set up. The middle (black) lines in 
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FIG. 7. (Color online) Cubic MnO: RIXS spectra for Mn 2+ 
with the 2p 6 3d 5 ground state. Full lines correspond to the 
vertical polarization of the incoming photon, dotted lines to 
the horizontal one. Top (blue) curves: raw results of the code. 
Middle (black) curves: experimental results Ghiringhelli et 
aZ.— . Lower (red) curves: results of the code with rescaling 
as given in the text, (a) and (b) are respectively taken at 0.7 
eV and 3.7 eV above the L3 edge. 



both panels are the experimental measurements^. The 
top (blue) lines are the raw results of the code, where a 
constant core-hole broadening of 0.2 eV and a constant 
final state broadening of 0.1 eV was taken. Although the 
raw spectra look somewhat shifted compared to the ex- 
periments, many of the features are already present. One 
notes in particular that the relative intensities for the two 
polarizations is reproduced. The bottom (red) curves in 
(a) and (b) are the spectra obtained after applying the 
scaling parameters S cou i=0.85, S x t a i=1.4 and S soc =l.l. 
The rescaled curves are in excellent agreement with the 
experiments and compare favorably with the fitting done 
with the group-theory crystal field approach in Ref. l39l . 
The latter fitting gave a value for the effective CXS of 
1 eV. In order to get an estimate of the strength of the 
crystal field in the present approach, the ground state 
multiplet calculation is repeated with only 1 d electron, 
in the non-relativistic limit. In that case the e g and ti g 
levels clearly define CXS. Without scaling the raw result 
gives CXS=0.74 eV. When applying the scaling param- 
eters CXS becomes 1.03 eV, an unsurprising result since 
the scaling parameters are deducted from the same fit as 
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FIG. 8. (Color online) Calculated MnO RIXS spectra for 
Mn 2+ with the 2p 6 3d 5 ground state; the same scaling param- 
eters as in Fig. [7] were applied. The incoming photon has an 
energy 3.6 eV above the L3 edge, and 6i n is in all cases along a 
Mn-O bond, (a) Results with averages around b ou t; top (red) 
line: cubic MnO with b out as in Ref. [H^; lower full (green) 
line: monoclinic MnO with b out nearly perpendicular to 6i n ; 
lower dashed (red) line: scaled cubic equivalent (see text), (b) 
Results with given e out ; top (red) lines for cubic, lower (green) 
lines for Mnl in monoclinic MnO; full lines: e ou t || ei n ; dotted 
lines: e out _L 6i n (only nearly perpendicular for monoclinic). 



in Ref. l39l . 

A high-resolution measurement in a slightly different 
geometry and taken at a different beamline was published 
three years later—. Fig. Uta) (top red line) shows the 
calculated spectra for this particular experimental set up, 
with the same scaling parameters as used for Fig. [7] The 
calculated spectra reproduces faithfully the features of 
the experimental measurement between -2 and -4 eV. 

As the present approach is very well suited to explor- 
ing the impact of different polarization directions on the 
spectra, we demonstrate in Fig. [5] (b) (top red curves) 
and Fig. 0(a) some of their effects in the cubic structure. 
To that purpose hypothetical spectra for selected incom- 
ing polarizations £i n and outgoing polarizations e out are 
drawn. 

The top (red) curves of Fig. [5] (b) show the calculated 
spectra of the cubic MnO for ej n along a Mn-0 bond; 
e ou t is either parallel (full line) or perpendicular (dotted 
line) to 6i n . As expected, the intensity of the elastic peak 
is weaker when ej n _L e out ; also in this case the peaks at 
around -3.6 eV are enhanced, whereas those at -2.2 and 
-4.1 are severely reduced. 

Fig. [9] (a) shows the effect of keeping e in _L e out , 
but with neither directions along any symmetry axis of 
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FIG. 9. (Color online) Calculated MnO RIXS spectra for 
Mn 2+ with the 2p 6 3d 5 ground state; the same scaling param- 
eters as in Fig. [7] were applied. The incoming photon has an 
energy 0.7 eV above the L3 edge, (a) Cubic MnO; 6i n _L e out , 
with neither ei n nor e out along a bond, (b) Top (red) curve: 
cubic MnO with e; n and b out in two arbitrary directions, and 
with the respective two directions swapped. Middle curves: 
Monoclinic MnO with various ei n and b ou t as given in the 
text (full green, maroon dashed, turquoise dash-dash-dotted 
lines), and with the directions of the first two swapped (green 
dotted, maroon dot-dot-dashed lines). Lower curves: differ- 
ences between curves above, as described in the text. 

the structure. The calculations for the full (turquoise) 
and dashed (maroon) lines were performed with ej n cho- 
sen deliberately away from any bonds, namely in the 
(0.4,0.5,0.77) direction, with respect to the unit cell vec- 
tors. The outgoing polarizations e out for those two lines 
are perpendicular to e- ln , being respectively in the (0.78,- 
0.62,0) and (0.35,0.690,-0.63) direction. For comparison, 
the dotted (red) curve shows ej n _L e out with ej n along a 
bond. In this situation any outgoing polarizations per- 
pendicular to it will produce the RIXS curve shown as 
the dotted line. When the incoming polarization is away 
from the symmetry axis, slightly different curves are pro- 
duced even though the outgoing polarization is at 90° 
from the incoming one. 



B. Low-temperature MnO 

MnO exhibits a phase transition^ from paramagnetic 
to antiferromagnetic state at 118 K. A neutron scatter- 
ing study^i of MnO at 10 K shows that the distortions 
accompanying the phase transition transform the struc- 
ture into a monoclinic one. Having established in the last 



section the sensitivity of RIXS to the different polariza- 
tion directions in the cubic case, the question arises as 
to what extent these deformations could be picked up by 
RIXS. 

Three different Mn sites were identified in the mono- 
clinic structure, all of them in distorted octahedral en- 
vironments with 6 oxygen atoms as neighbors. The 
three distorted octahedra, although not exactly aligned, 
have similar orientations. In the cubic structure, drastic 
changes in the RIXS curves are observed for the outgo- 
ing polarization parallel and perpendicular to that of the 
incoming light, Fig. [Sfb), top (red) lines. Effects of a 
similar magnitude are expected to occur in the mono- 
clinic structure. 

First the hypothetical spectra with selected incoming 
polarization ej n and outgoing polarization e out are con- 
sidered, Fig. \8\ (b). The lower (green) curves of Fig. [8] 
(b) show the calculated spectra for Mnl only, with ej n 
along the Mnl-Ol bond; the full line has been obtained 
for e out || e in , the dotted line for e out along the Mnl- 
02 bond, which is nearly perpendicular to the Mnl-Ol 
direction. The comparison between the top (cubic) and 
lower (monoclinic) curves of Fig. [Sfb) shows that the dis- 
tortions are causing only minute changes to the spectra. 
One notes that the peak originally at -2.2 eV is shifted 
to the right by 80 meV, and that at 4.1 to the left by 50 
meV, accompanied by some minor changes in height. 

The predicted spectrum for monoclinic MnO for a 
given 6i n and a detection direction 6 out is plotted in Fig. 
[5Ja) as the lower full (green) curve. For an optimal com- 
parison with the cubic case, ej n has been taken along the 
Mnl-Ol bond, and 6 ou t along the nearly perpendicular 
Mnl-02 direction. The curve is the result of an aver- 
age over the three Mn sites. In order to show the effect 
of distortions, the spectrum for the cubic structure has 
been recalculated with the point charges distances uni- 
formly scaled down (dashed red line), so that the crystal 
field of the cubic structure matches the average of the 
distorted one. In effect, the cubic Mn-0 distance d has 
been changed so that d~ 3 is the same as the average 
(r~ 3 ) of the distorted octahedra, so that Mn-0 distance 
is 2.212 instead of 2.22 A. The differences between the 
full (green) line and the (red) dashed line therefore ex- 
emplify the expected differences between the cubic and 
distorted structures. The same types of differences be- 
tween the cubic and monoclinic can be noted as above in 
the discussion of Fig. [H^b). 

A more stringent test for the distortions can be worked 
out by noting that in the cubic structure, the same RIXS 
spectrum is obtained by exchanging the direction of the 
incoming polarization Cj n with that of the unpolarizcd 
detection direction b out , for any chosen directions. An 
example is plotted in Fig. |9](b), top (red) curve. Two di- 
rections were picked randomly, di = (0.18, —0.40, —0.90) 
and d 2 = (0.56, —0.31, —0.77) (with respect to the unit 
cell vectors). The same curve is produced by taking 
ein=di with b out =d 2 , and e in =d 2 with 6 out =di. This 
is not the case for the distorted structure, for which 
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FIG. 10. (Color online) RIXS spectra for Ni 2+ with the 
2p 6 3d 8 ground state in NiO. Top curve (blue): raw results 
of the code. Middle curve (black): Experimental result of 
Ghiringhelli et a/.—. Bottom curve (red): result of the code 
after scaling (see text). 
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FIG. 11. Low-lying multiplets for Dy 3+ in the 3d 10 4/ 9 ground 
state, (a) and (b) Calculated levels in zero crystal field with 
respectively S cou i=S SO c = l, and S cou i=0.85, S SO c=0.95, and as- 
sociated multiplicity, (c) Ion spectroscopy experiments^ and 
associated J value, (d) Calculated multiplets in crystal field 
(see text) with S CO ui=0.85 and S SO c=0.95. 



slightly different curves are calculated when the direc- 
tions arc swapped. Some examples are shown as the 
middle curves in Fig. [9](b). The full and dotted (green) 
curves are produced with di along the Mnl-Ol bond, 
and d2 along the monoclinic crystallographic axis a; the 
full curve has ej n =di with 6 out =d2, the dotted curve 
£in=d2 with 6 out =di. The difference between these two 
curve is plotted as full (green) line in the lower part of 
the figure. A greater difference is observed by taking di 
along a, and d2 perpendicular to the crystallographic di- 
rections a and 6, which we call d . The relevant curves are 
drawn as (maroon) dashed and dot-dot-dashed, with the 
difference plotted as dashed line. The middle (turquoise) 
dash-dash-dot curve was calculated with ej n along 6, and 
6 out along a. The difference between this spectrum and 
one with e; n along d and b out along a is shown in the 
lower part of the figure as the (turquoise) dotted line. 
Here the difference is greater in the lower part of the 
spectra. It is expected that those tiny differences, in this 
case smaller as 10% of the intensity of the curve, will be 
very difficult to detect experimentally. 



C. Cubic NiO 

Another example where a simple crystal field model 
can be applied is at the Ni L-edge of NiO^. The high- 
resolution experiment of Ghiringhelli et al£Q is repro- 
duced as the middle (black) curve in Fig. [TUJ Adopting 
the set up described in Ref. the calculations were car- 
ried assuming cubic symmetry for NiO with the shortest 
Ni-0 distance at 2.09 A. Ni in NiO as the valance Ni 2+ 
and the ground state is therefore 2p 6 3d 8 . The top (blue) 
line in Fig. \TU\is the raw result of the code. The bottom 
line (red) is obtained by applying the scaling parameters 
S C oui=0.7, S xt ai=2.0 and S SOC =0.915. A constant core- 



hole broadening of 0.1 eV and final state broadening was 
applied to both calculated curves. Although again the 
raw curve contains the right features with the correct hi- 
erarchy of peaks, a fairly strong rescaling of the crystal 
field strength had to be applied in that case to bring the 
curve in good agreement with the experiment. The CXS 
value can be obtained in a similar fashion as for MnO by 
placing a single d electron in the 3d shell and perform- 
ing the calculation in the non-relativistic limit. The raw 
multiplct result gives an unsealed CXS of 0.51 eV and a 
scaled CXS of 1.03 eV. The latter value is in complete 
agreement with the group-theory approach of Ref. l52l . 



D. Dy phthalocyanine 

Lastly, we use our code to predict the crystal-field ef- 
fects on the RIXS spectrum of a lanthanide compound, 
the "double-decker" phthalocyanine complex^, with for- 
mula C64H32DyiNi6. The multiplet levels for this struc- 
ture have been determined by Ishikawa et al by com- 
bining NMR measurements and ligand-field parameters 
modeling^, thereby providing a point of comparison with 
the ground state multiplets calculated by our code. In 
the selected structure a trivalcnt Dy atom sits between 
two planar ligands composed of the N, C and H atoms. 
Dy 3+ is in the ground state 3d 10 4/ 9 . 

Starting without a crystal field, the energy levels for 
Dy 3+ can be calculated and compared with those re- 
ported in the literature for Dy IV—. The latter are 
plotted in column (c) of Fig. [TTJ with their correspond- 
ing J value. The columns (a) and (b) show the lowest 
levels and their multiplicity, respectively with no scal- 
ing (a) and by reducing the intra-coulomb interaction by 
S CO ui=0-85 and the spin-orbit coupling by S SOC =0.95 (b). 
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The raw result of the code is in excellent agreement with 
the measurements (c) for the first levels, with a small 
shift upwards for the higher plotted values. The match 
is further improved by the application of rescaling. Two 
of the calculated levels (dashed lines in the figure) were 
not reported in Ref. HH; the other levels are in one to one 
correspondence with the measurements. 

As starting point for the crystal field charges, the oxi- 
dation states N -3 , C +4 and H +1 are a poor choice in view 
of the non-ionic character of this compound. A DFT cal- 
culation indicated that the charges are close to being a 
tenth of the oxydation state values^. As an alternative 
to scaling the charges down, one can set the crystal field 
scaling parameter S x tai to 0.1. Column (d) show the first 
54 levels calculated with this crystal field and the scaling 
parameters S CO ui=0.85 and S SOC =0.95. One notes that 
the 16-fold degenerate ground state has split under the 
influence of the crystal field in a series of low-lying, two- 
fold degenerate states, all very close to zero. Ishikawa et 
al find the first level at 0.004 eV and the last at 0.068 eV; 
the present code is in close agreement with respectively 
0.008 and 0.059 eV. 

The XAS and RIXS spectra for the Dy phthalocyanine 
complex calculated with these scaling parameters are 
plotted in Fig. [TJ The XAS in Fig. [TJ shows the 
M 5 part of the spectra and indicates the energies at 
which the RIXS spectra were calculated. The following 
incoming polarization directions were used: along 
the crystallographic axis a (full red line), along c 
(dash-dotted violet line), and along the C4 axis (green 
dot-dot-dashed line) of the molecule. The (black) dotted 
line shows the result in absence of crystal field. The 
energies B and C are respectively at 2.3 and 4.9 eV 
from the point A. The RIXS spectra were calculated 
for various incoming polarizations e; n and detection 
directions b ou t- The full (red) lines were generated for 
£in = b out both along the crystallographic axis a. The 
dashed (turquoise) curves also had e; n along a, but b out 
along c, while the dash-dotted (violet) lines took €i n 
along c, and 6 ou t along a. Setting e- m = b out along 
the C4 axis of the molecule gives the dash-dot-dotted 
(green) curves. The RIXS intensity with no crystal 
field (black dots) were scaled down by a factor 8 for an 
optimal comparison with the other curves. This factor 
reflects the respective degeneracy of the ground state 
levels with and without the application of the crystal 
field. The lowest lying states are absorbed by the elastic 
peaks, but the states between 0.4 and 1 eV are visible. 
Another feature is the predominance of the peak at 2.3 
eV at energy B, then of the peak at 2.7 eV at energy 
C, for the incoming polarization along a. Although it 
is likely that the chosen scaling parameters would need 
some adjusting, testing with parameters indicated that 
the major effect on the RIXS spectra of Fig. [T^] will 
appear to produce a shift of the peaks between 2 and 4 
eV along the energy axis, rather than drastic changes in 
polarizations behavior. 




Energy loss (eV) 



FIG. 12. (Color online) Calculated XAS and RIXS spectra 
for Dy 3+ with the 3d 10 4/ 9 ground state in Dy phthalocya- 
nine. XAS: with polarization along a (full red line), c (dash- 
dotted violet line), and C4 axis (dash-dot-dotted green line); 
(black) dotted line are calculated in the absence of crystal 
field. RIXS: full (red) lines: ei n = b out || a, dashes (turquoise): 
6i n || a, b out II c, dash-dots (violet): 6i n || c, b out || a, dash- 
dot-dots (green): ei n = b out along C4 axis of the molecule; 
the dotted (black) lines have been produced with e- m = bout 
in the absence of crystal field, and have been scaled down by 
a factor 8. 

In this section RIXS calculations on cubic MnO and 
NiO were tested against measurements. The raw curves 
exhibited the correct relative behavior, and very good 
agreement with experiments could be achieved by apply- 
ing scaling parameters. The level of changes in the spec- 
tra caused by distortions could be evaluated for the low 
temperature MnO . Finally, the code was applied on an / 
electron compound, Dy phthalocyanine. After establish- 
ing by comparisons with measurements that reasonable 
multiplets levels are obtained for the ground state, RIXS 
predictions were made. 



VI. CONCLUSION AND OUTLOOK 

In summary, we have described an easily accessible 
approach for the computation of multiplet spectra aris- 
ing from soft X-ray spectroscopies of narrow band solids. 
The required input data is 

• The element being probed (for example Ti) 

• The core and valence configuration of the ground 
state for the shells involved in the transitions (for 
example 2p 6 id 2 ) 
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• The crystal field as a list of positions and charges 

The charges can be viewed as semi-empirical, with for- 
mal charges or oxidation states as a starting guess. The 
central field part of the present theory is done with a 
Dirac relativistic implementation for a simple density 
functional, and yields the spin-orbit splitting from first 
principles. The electron-electron interaction is calculated 
explicitly for the electrons of the open shells. The crystal 
field is obtained as the resulting electrostatic field gener- 
ated by the surrounding point charges. Any crystal field 
of arbitrary and unspecified symmetry can therefore be 
entered without expert knowledge of group theory no- 
tation. Moreover, this approach allows the polarization 
directions to be given in a straightforward manner with 
respect to the positions of the crystal field ions. Within 
the limitations of the model and the scope of the approxi- 
mations made, the method constitutes a "first principles" 
approach for calculating multiplets states as well as XAS 
and RIXS spectra. In the latter cases the core-hole life- 
time is also required as an input . 

Our examples show that the raw, "first-principle" re- 
sults provide a good estimation of the multiplet structure; 
the correct experimental trends are reproduced in those 
experiments or part of the spectrum where the crystal 
field plays a dominant role. The model of the single atom 
in an ionic environment limits however the agreement 
with experiments. This can be corrected to some extent 
by the introduction of semi-empirical scaling parameters 
that account for screening effects, hybridization and cou- 
pling to energy bands that are not originally within the 
scope of the model. The scaling parameter for the intra- 



atomic Coulomb interaction varies typically between 0.7 
and 1; a 10% adjustment of the spin-orbit coupling has 
sometimes been necessary; in the examples discussed in 
this paper the value of crystal field scaling parameter has 
ranged from 1.0 and 2.0. These corrections to the raw 
result reflect the fact that in a bonded solid, the radial 
functions extend further than anticipated by the neutral 
atom calculation. The crystal field correction may also 
include a global adjustment to the point charge values. 

Although the use of scaling parameters leads to ex- 
cellent fits to experiments, it is necessary to build up 
the model and include non-local contributions into the 
spectra. The crystal field model for any symmetry con- 
stitutes the first step towards a more complete approach 
to multiplets for core spectroscopies. 
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